Phase Equilibrium of Binary Mixtures in Mixed Dimensions 
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We study the stability of a Bose-Fermi system loaded into an array of coupled one-dimensional (ID) "tubes", 
where bosons and fermions experience different dimensions: Bosons are heavy and strongly localized in the ID 
tubes, whereas fermions are light and can hop between the tubes. Using the ^^'*Yb-''Li system as a reference, 
we obtain the equilibrium phase diagram. We find that, for both attractive and repulsive interspecies interaction, 
the exact treatment of ID bosons via the Bethe ansatz implies that the transitions between pure fermion and any 
phase with a finite density of bosons can only be first order and never continuous, resulting in phase separation in 
density space. In contrast, the order of the transition between the pure boson and the mixed phase can either be 
second or first order depending on whether fermions are allowed to hop between the tubes or they also are strictly 
confined in ID. We discuss the implications of our findings for current experiments on ^'^Yb-®Li mixtures as 
well as Fermi-Fermi mixtures of light and heavy atoms in a mixed dimensional optical lattice system. 
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I. INTRODUCTION 

The quest for lower temperatures in ultracold gases has lead 
to the development of many ingenious techniques to cool sev- 
eral atomic and molecular species, hi particular, the explo- 
sion of activity concerning ultracold Fermi gases has become 
possible largely owing to the success of sympathetic cool- 
ing, which allows to efficiently cool fermions by mixing them 
with bosons |[T|. At the same time, this procedure has stim- 
ulated the investigation, both experimental and theoretical, of 
Bose-Fermi mixtures. Here, the possibility of tuning the inter- 
species interaction strength using Feshbach resonances lEJ^I, 
has led to the exploration of many interesting phenomena such 
as collapse and phase separation |5, 6|, as well as boson- 
mediated Cooper pairing |1|. Furthermore, Feshbach reso- 
nances can also be used to generate heteronuclear molecules, 
which can exhibit large electric dipole moments. This opens 
the interesting possibility of studying dipole-dipole interac- 
tions in quantum degenerate gases [7-:- .13] . 

Meanwhile, the advent of optical-lattice confinement ID 
[141 has turned ultracold atomic gases into unique environ- 
ments where to simulate and understand strongly correlated 
phenomena relevant to condensed matter systems. This has 
been made it possible, for example, by confining the atomic 
clouds in reduced dimensions, such as a one dimensional (ID) 
array of two-dimensional (2D) planes or a 2D array of ID 
tubes. Whereas the former has enabled the study of interest- 
ing phenomena occurring in two-dimensions lfT5l[T6ll . the lat- 
ter has provided us with an amazingly tunable tool to explore 
the physics of interacting ID quantum systems lfT7H20l . of 
which it is much more difficult to find faithful realizations in 
a more conventional condensed matter context. 

Optical lattice confinement has also allowed to envisage the 
realization of new types of quantum systems. One such exam- 
ple, analyzed in this work, is provided by mixtures of inter- 
acting particles in mixed-dimensional lattices. In recent years, 
these systems have attracted an increasing amount of theoret- 
ical attention Ii21ii22j . and very recently they have been also 



experimentally realized 11231 124 1. Besides its intrinsic inter- 
est as a new category of quantum many-body systems, they 
may also offer advantages for reducing few-body losses and 
enhancing stability in strongly interacting regimes f25\. 

In recent years, there has been a growing interest in un- 
derstanding the properties of mixtures of ultracold Bose 
and Fermi gases. In particular, binary mixtures of bosons 
and spin-polarized fermions have been studied in three- 
dimensions (3D) IH IS Eim, 2D EO] and ID EBtia. 
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FIG. 1: (Color online) Schematic representation of the system stud- 
ied in this work, namely a mixed dimensionality Bose-Fermi system. 
The Bose-Fermi cloud is loaded in an anisotropic optical lattice that 
produces a two-dimensional A^^ x N^ array of one-dimensional (ID) 
tubes of length L. The fermions are light enough to hop between the 
tubes whereas bosons, which are assumed to be much heavier, are 
confined in the ID tubes. The bosonic and fermionic Wannier func- 
tions {ip-R.(r±) and TOR(rx), respectively) are also displayed. 



Note that the ID geometry has a special relevance because 
of the central role played by quantum fluctuations and the fact 
that there is neither broken continuous symmetry nor, conse- 
quently, off-diagonal long-range order. The equilibrium phase 
diagram of ID Bose-Fermi mixtures has been considered by 
many authors ll3T1435l . while the case of a Bose-Fermi system 
in an anisotropic optical lattices was studied by one of us in 

Ref. ma. 

In addition to dimensionality issues, the significant interest 
in studying the stability of binary mixtures comes also from 
the possibility of tuning the Bose-Fermi scattering length via 
the Feshbach resonance mechanism. Indeed, in the repulsive 
interaction regime, the spatial overlap between bosons and 
spin-polarized fermions is reduced, thus ensuring the stability 
of the system |36, 37 1. When the repulsion is increased, the 
two components tend instead to phase separate, rather than 
uniformly mix: In the particular case of a 3D geometry, phase 
separation occurs either between a mixed phase and a purely 
fermionic phase or between two pure phases [6]. In the regime 
where the interaction between bosons and fermions is attrac- 
tive, a significant reduction of the interatomic distance can 
lead to a collapse of the mixture, because of three-body re- 
combination processes ll38l . As discussed later, the stability 
of the mixture towards a collapsed phase can be enhanced in 
one-dimensional geometries |25|. 

In this article, we study the phase stability of a Bose-Fermi 
mixture embedded in a mixed-dimensional optical lattice of 
an array of one-dimensional tubes |39, 40| (see Fig.flll. The 
mixed-dimensionality comes from the fact that, while bosons 
are longitudinally confined along the tubes and strictly move 
in ID, fermions are not constrained to ID and are allowed 
to hop between tubes in the transverse directions. This as- 
sumption can be justified based on the fact that many real- 
izations of Bose-Fermi mixtures consider bosonic species that 
are much heavier than the fermionic ones. We do also assume 
that, while bosons interact with each other, as well as with 
fermions, the fermionic component is polarized in a single 
hyperfine state and thus is non-interacting at very low tem- 
peratures. 

As already mentioned, one motivation to study this ge- 
ometry is that, by confining bosons in ID, the mixture sta- 
bility is enhanced ll25l . Note that, even for purely bosonic 
gases, the spatial overlap between bosons has been experi- 
mentally shown to be strongly suppressed by the strong corre- 
lations emerging in ID 141]. In addition, mixed-dimensional 
systems allow to study interesting few-body |42| and many- 
body II22I l43l phenomena. We should stress that the re- 
sults reported here are relevant to several experimental re- 
alizations. The simplest is, for example, a mixture of light 
fermions and heavy bosons. In particular, in order to make 
contact with ongoing and future experiments ll44H46ll . we ex- 
plicitly consider a mixture of ^Li (a light fermion) and ^^'^Yb 
(a heavy boson) atoms. An alternative realization could be 
given by an originally imbalanced mixture of fermions in two 
hyperfine states |47|: Bosons are then formed by associating 
fermions into Feshbach molecules |47|, leaving out from the 
pairing the spin polarized excess majority fermions. When 
the fermions belong to the same atomic species, the Fesh- 



bach molecules have twice the mass and twice the polarizabil- 
ity of the fermionic atoms. Therefore, it should be relatively 
easy to confine the bosons in ID by loading them into a two- 
dimensional optical lattice. On the same lattice, the remaining 
majority fermions behave as the light component. 

We find that, since bosons are confined to strictly ID, they 
can undergo fermionization. This means that, as they mix with 
fermions, they form a Tonks-Girardeau gas whose energy per 
unit length grows as third power of the lineal density. As a 
result, we find that the nature of the transition from the pure 
Fermi gas to a mixture is always first order, implying phase 
separation in density space. Note that this is a very different 
result from the one obtained by assuming that bosons form 
a (quasi-)condensate, with the energy density growing as the 
square of the boson lineal density. The latter situation applies 
either to a high density Bose gas in ID or to a gas of bosons 
that can hop in 3D. We also find that, in the mixed dimen- 
sionality lattice, the transition between pure boson and mixed 
phases is continuous, while it becomes first order when the 
fermions are also confined in ID. 

The rest of this article is organized as follows: In Sec.[ll|we 
introduce the model for the Bose-Fermi mixture and discuss 
methods and approximations employed to derive the system 
phase diagram. In particular, in Sec. II A we explain how 
a mean-field approximation is applied solely to the boson- 
fermion interaction, whereas the boson-boson interaction is 



treated non-perturbatively using the Bethe ansatz in Sec. II B 
In Sec 



III 



we derive the phase diagram and interpret the ori- 
gin of the discontinuous character of the transitions between 
the pure fermion and mixed phases in terms of an expansion 
of the free energy for small values of the boson density. We 
first describe the results obtained for the case of mixed dimen- 



sions in Sec. Ill A and later for the pure ID limit in Sec. Ill B 



Finally, in Sec. IV we present the main conclusions of our 
study and discuss the limitations of our approach, as well as 
some directions for future work. Some technical aspects of 
our derivations are provided in the appendices. 



II. MODEL 

We consider a mixture of interacting bosonic (B) and 
single-component fermionic (F) atoms described by the fol- 
lowing Hamiltonian (in ^ = 1 units): 



H — Hb + Hp + Hbf 

V2 
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where the density operators are ^^(r) = ^Jj(r)^^(r), with 
a = B,F. We have approximated all interaction potentials 
by contact interactions, which are parameterized by an s-wave 



scattering length Uap: 
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being ttibb = ttt-b and tubf = 2mBmp/{mB + mp). For 
thermodynamic stabiHty reasons, the bosons are assumed to 
repel each other (i.e., ubb > 0). The value of ubb can 
be tuned by, e.g., controlling the strength of the transverse 
confinement |48|. The sign of the Bose-Fermi interaction 
strength is determined by ubf^ which can be controlled inde- 
pendently from asB by resorting to an inter-species Feshbach 
resonance. Below we consider both the repulsive (ubf > 0) 
as well as the attractive (gbf < 0) case. At ultracold temper- 
atures, interactions between identical fermions can be safely 
neglected. 

The Bose-Fermi mixture is loaded into an anisotropic op- 
tical lattice formed by a 2D Ny x N^ square array of ID 
"tubes" of length L directed along the cc-direction and equally 
spaced by a distance d (see Fig. [Til. This can be described by 
an optical potential of the form VQ(rj^) = Vf^[sm^ (ny/d) + 
siii^ {nz/d)], where rj^ = {y,z). The strength of the optical 
potential VJ,^ (^o^) i^ determined by both the laser intensity 
and the boson (fermion) atomic polarizability, allowing the 
possibility of mixed dimensionality for the mixture. In par- 
ticular, we assume that the bosons are tightly confined in ID 
tubes and thus move strictly in ID, while fermions can hop be- 
tween the tubes. We will derive the thermodynamic phase dia- 
gram for this geometry, thus neglecting the harmonic confine- 
ment. By making use of the local density approximation, in- 
formation about the experimentally relevant trapped case can 
be extracted from the homogeneous phase diagram plotted in 
chemical potential space. 

Because the bosons in the mixture are assumed to be 
more massive than fermions and/or to experience a more 
confining lattice potential yB(rj^), they are tightly confined 
along the "tubes", a configuration often refereed to as a two- 
dimensional optical lattice d [18] [20l gO] EH [SO) . Thus, the 
field operator *I'^(r) can be expressed in terms of Wannier 
functions (/9j^(rj^) localized at the tube site R = (iy,i^)d 
(see, e.g., ID and references therein) and the tube boson oper- 
ator ^^rW: 



R 



(3) 



The Wannier functions form an orthonormal basis. By ne- 
glecting the interactions between tubes with R 7^ R', we can 
rewrite the bosonic Hamiltonian Hb in ([T]i as a sum of decou- 
pled ID Hamiltonians, 
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PBuix) 



*br(2;), (4) 



where pBuix) = ^ b'r(^)'^ Biii^) ts the single-tube boson 
density operator, g^^ is the one-dimensional boson coupling. 



which, for weak boson-boson interaction, takes the form: 

9bb = 9BB / dr±\ip^{r±)\^ , (5) 



and /ijj is the ID boson chemical potential: 



pW = / drj_ip^{r±) 



2mB 



-VB{r±_) + fiB 



<y5R(r-L) ■ 



Note that, as we will see later, for strong boson-boson interac- 
tions, the expression of the ID boson coupling g^g is instead 
given by Eq. (|28]l |48 1 rather than Eq. ^. 

In contrast, we assume the fermions to be more weakly con- 
fined than bosons along each tube due to their smaller mass 
and/or a weaker optical potential. Yet, the lattice confinement 
is strong enough so that the description of the Fermi field in 
terms on the lowest Bloch band t/ik^ (r_L) is accurate and we 
can expand: 



*f(i-) 
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y^'gifex 



<_l(i"±)/i 



k ' 



(6) 



where k = (fc^,kj^). Here, whereas the motion along the 
X direction is free, the motion in the transverse directions 
rj^ — {y, z) is described by (^k^ (rj_), which is a Bloch wave- 
function belonging to the lowest Bloch band characterized by 
a crystal momentum kj^. Projecting the fermion Hamiltonian 
onto this band yields: 



Hf^J2 [^(k) - Pf] Ilk , 



where the Fermion dispersion reads as: 



e(k) 



e(kj 



2mF 
e(kj_) = 2i [2 — cos{kyd) — cos{kzd)] 



(7) 

(8) 
(9) 



Thus, summarizing, in the geometry studied here, the 
bosons are tightly confined to move in ID, whereas the 
fermions can hop between the tubes, although the optical lat- 
tice potential does affect thek dispersion relation. Under these 
conditions, it is known ll20l that, at low temperatures, the 
bosonic atoms loose their individuality and the low-energy 
long-wavelength excitations are ID phonons. For arbitrary 
values of g]^ > 0, the ground state properties of such an 
interacting bosonic gas are exactly described by the Bethe 
ansatz solution obtained by Lieb and Liniger iBTI . The lack of 
individuality and the highly correlated behavior brought about 
by the ID confinement calls for a treatment of the problem 
that treats the boson-boson interactions in a non-perturbative 
way. Since a mean-field approximation fails to capture the 
fundamental bosonic correlations in ID, we apply it only to 
the interactions between the fermions and the bosons, as we 
explain in what follows. 



A. Bose-Fermi interaction: mean-field 

In order to render the above model tractable, we apply a 
mean-field approximation to the Bose-Fermi interaction term. 



We do this in such a way that the different ID tubes become 
decoupled at the expense of introducing self-consistent shifts 
of both the boson and fermion chemical potentials. To this 
end, we first observe that the tight confinement of the bosons 
in ID allows us to neglect the overlap between Wannier func- 
tions localized at different ID tubes (see Fig. [TJ and thus we 
can approximate the density operator of the bosons as 

Ps(r) - Yl \VRi-^±)\^PBRix) . 

R 



it is obtained in the limit where the transverse confinement for 
the bosons is tight (see App. |A] for the details of the deriva- 
tion). Also, we have introduced the following lineal densities: 



P° = 



NL 



(15) 



where N = NyN, 



is the total number of ID tubes. Thus, 
within this mean-field approximation, the system Hamilto- 
nian ([T]l can be written as 



In this limit, the Bose-Fermi interaction term Hbf of the 
Hamiltonian ([TJ can be written as 



Hbf ^ gsF ^ / dxpBR{x)pFR{x) 

R -^ 



(10) 



where pfr{x) is a projection of the 3D Fermi density operator 
/5i?(r) on the R-th tube: 



Pfr{x) = / dr_L|v5R(r_L)|Vi=^(r) 



(11) 



We note that this approximation does not decouple the dif- 
ferent tubes yet because, even if bosons cannot hop from one 
tube to another, the interaction between bosons belonging to 
different tubes is mediated by the hopping fermions, i.e., the 
operator pFRix). However, if we rely on a mean-field ap- 
proximation to replace the operator pfr{x) by its expecta- 
tion value (which, as shown in detail in App. [A} is a constant), 
the different tubes become decoupled, which allows for a so- 
lution of the model introduced above. We emphasize again 
that this kind of mean-field approximation is different from 
the standard treatment (for example employed in Ref. |25|), 
where also the boson density in the boson-boson interaction 
term is replaced with its expectation value. Instead, here the 
boson interaction is treated non-perturbatively using the Bethe 
ansatz, emphasizing that the fundamental entities subject to 
the mean-field interaction are the ID tubes and not the bosons 
themselves. 

Thus, in Eq. ( [TO] l, we write the density operators PaR{x) 
as their quantum averages plus fluctuations, i.e. PaR{x) = 
{PaR{x)) + SpaR{x). Hcncc, the mean-field approximation 
is obtained by substituting these expressions into Eq. ( [TO] i and 
by neglecting the second order terms in the fluctuations, which 
leads to: 



H^ 



BF 



H 



mf _ 



X^ 



dx [-{pBRix)) {pfr{x)) 



BF - 9BF ^_^ 
R 

+Pbr{x){pfr{x)) + {pbr{x)) Pfr{x)] ■ (12) 

In absence of harmonic confinement along the tubes, transla- 
tional invariance along the x-direction requires that the aver- 
ages 



(PBRix)) 



P% 



{pfr{x)) ^ A p^p (13) 



are constants independent on the tube index R. Here, the con- 
stant 

Ac,Njdr^\^^{rj_)\^\<PuAr±)W (14) 
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H ~ H"'^ = H 
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E 
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gBFAp%p^p , (16) 



where Hr, is defined as the bosonic Hamiltonian from 



's 



Eq. Q with the chemical potential shifted as p^^P 
cjBFAp^p: 
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and -ff™ is the fermion Hamiltonian from Eq. ^ with 
chemical potential shifted as pp —>■ pp — gBFAp''^: 



H-^ =. J2 [e(k) -PF + 9bfAp%] flf^ 



jmf _ 



(18) 



We would like to stress that we are not applying a mean-field 
approximation to the boson-boson interaction term p^j^(x). 
On the contrary, as shown in the next section, we shall treat 
this term exactly using the Bethe ansatz solution of Eq. ( [T7| ) 
due to Lieb and Liniger Ii5 1 J . 



B. Zero-temperature free energy 

Starting from the mean-field Hamiltonian defined by equa- 
tions ([T6|, ( [T7] i, and ( fTSJ ), we can evaluate the grand-canonical 
free energy density at zero temperature. Note that, by virtue 
of the mean-field approximation and the translational invari- 
ance, the bosonic '^ ^-^^ and fermionic /j. field operators in 
Eq. ( [T6] l have become decoupled and therefore we can sepa- 
rate the bosonic and a fermionic contributions to the free en- 
ergy potential / = f{p^, pp, p^), which can be written as: 



/ 



ip%f 
2mB 



e(7) 



ID 

Ms Pb 



1 

N 



E 



l3 

Up. 



Jk^) 



Sirmp 



(19) 



where BZ stands for the 2D Brillouin zone, i.e., the region 
of kj^ -space where \ky^z\ < n/d. In evaluating this expres- 
sion, we have relied upon the Bethe ansatz solution 12011511 of 
the interacting ID boson Hamiltonian, H^g ^T7) . It is worth 
noting that the constant term gsFAp'^p'^ in Eq. ([T6| cancels 
exactly the mean-field shift of the bosonic chemical potential, 
P]F -^ P]F - 9bfAp1„ in Eq. ^. 



In Eq. ( [T9] l, the dimensionless function 6(7), where 7 — 
rUBg^/^B, is determined by numerically solving the fol- 
lowing system of coupled integral equations: 



e(7) 
2TTg{u) 






1 + 21 



duu g{u) 



du' — 

-1 [u 



g{u') 



u'Y + P 



(20) 
(21) 



where I — ^ j_^du g{u). The fermionic contribution to the 
free energy can be expressed as an integral over the transverse 
momentum, which leads to the last term in Eq. ([T9|, where we 
have defined: 



fc_Fx(k_L) = Rey 2to_f [^_f - gspAp^ - e(k_L)] . 

Finally, the thermodynamic grand-canonical free energy 
density is obtained by finding the global minimum of the po- 
tential f{fi^,fiF, Pb) with respect to the boson density p^: 



^il^B ,tJ-F)= min fipg ,Pf,Pb) 



(22) 



where ub denotes the equilibrium lineal boson density. Note 
that, since Bose-Einstein condensation is not allowed in ID 
interacting boson systems 1201 . the boson density cannot be 
regarded as the condensate density, i.e., the square of the 
condensate order parameter In addition, the equilibrium ID 
fermionic density can be evaluated from 



np 
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'n 



E 



kpxi^A 



(23) 



As explained in the next section, we can thus now evaluate the 
system equilibrium phase diagram either in chemical potential 
or in density space. 
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FIG. 2: (Color online) Zero temperature phase diagrams for a re- 
pulsive (aBF > 0, top panel) and attractive ubf < (bottom 
panel) Bose-Fermi mixture with fixed dimensionless hopping ampli- 
tude t — 37.8 and fixed interaction parameter ^ — 0.23 (the param- 
eters of this calculation correspond to a ^'^^ Yb-^Li system, see main 
text). The diagrams are plotted vs. the boson p}^ /Etb and fermion 
Pf /Etb chemical potentials, where Etb is the boson recoil energy 
— note that p^^ / Etb ^ 1, thus ensuring that the bosons remain 
confined to ID throughout. The thin solid black lines correspond to 
2nd order (i.e. continuous) phase transitions between either the vac- 
uum and the pure boson/fermion phases or between the pure boson 
and the homogeneous mixed phase. The thick solid (red) line cor- 
responds to a first order transition between the pure fermion and the 
mixed phase. Second and first order lines meet at at the tricritical 
point at Pf = Pb ~ (filled blue circle). The region with finite 
boson density is the light gray shaded region. In both cases phase 
separation can only occur between a mixed and pure fermion phases. 



III. PHASE DIAGRAM 

In this section, we obtain the phase diagram of the mixed- 
dimensionality system by minimizing the free energy intro- 
duced above in Eq. ^T9\ with respect to the boson density pg 
for fixed p^ and pp- The free energy f{p^J-', pf, p'b) also 
depends on several other parameters, such as gBF, gMs' ^'^'^ 
t, as well as the particle masses uib, rriF- Thus, it is conve- 
nient to simplify the description of the system by considering 
the following minimal set of four independent dimensionless 
parameters: 
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ID 
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2mBgBB 



(24) 
(25) 



where c = 2Am,p m^ \gBF\- Thus, the dimensionless in- 
teraction parameter 7 can we rewritten as 7 = (2^/5^)^^, 



with the dimensionless boson density given by p% = p%/c. 
In terms of the above dimensionless quantities, the free energy 
takes the form: 



i=iP%M,)-A-Pi-l, E %Sr. 



(26) 



k^esz 



where kFx{^i^ = Re^/2_F - sign{gBF)p% - ^O^-l) and 
e(kj_) = 2i[2 — cos{kyd) — cos{kzd)]. 



A. Mixed dimensions 

We explicitly consider here the case of mixed dimensions, 
while later in Sec. 



IIIB 



we will derive the phase diagram 
for the case of pure ID. We have numerically minimized / 
by fixing the value of the dimensionless interaction parame- 
ter, (, and hopping amplitude, t. In order to make contact 
with ongoing as well as future experiments |44-46J, we con- 
sider the specific case of a Bose-Fermi system consisting of 



a light fermionic species such as ^Li and a heavy bosonic 
species Uke ^^''Yb. When this system is loaded in a suffi- 
ciently deep 2D optical lattice, the large boson to fermion 
mass ratio {mB/rap — 29) is enough to suppress the hop- 
ping between tubes of bosons, while allowing fermions to hop 
between the tubes. This makes it possible to realize our initial 
assumption of mixed dimensionality for the system. Indeed, 
in the limit of a deep lattice, the fermion hopping amplitude 
in the tight-binding approximation of Eq. (|7ji can be expressed 
in terms of the optical potential strength V^ and the Fermi re- 
coil energy Erp, where Er^ = 2Tr'^ /{niaX'^) and X = 2d is 
the wavelength of the laser generating the optical lattice po- 
tential lU: 



4Erf 
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^-2y'V,f/ErF 



(27) 



For a laser wavelength A = 1064 run, the deep lattice con- 
dition is achieved by making Vq^ ~ AOErg (for this sys- 
tem V(f ~ 2VfP ISl EU). Furthermore, the boson-boson 
scattering (|2]) length has been experimentally estimated to be 
dBB = 104.9 qq ll53l (where ao is the Bohr radius). The ID 
interaction strength g^ can be obtained from |48 1: 
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where C ~ 1.0326 and is — 67.34 nm {6T\. For these pa- 
rameters, £f — 131.40 nm, and therefore, using Eq. ( |A5| l, 
A= 1.46017xl0i3m-2 Finally, we set [as f| = 13 an 1441. 
and allow for both positive and negative signs for asp, i-C- 
for the Bose-Fermi interactions to be repulsive or attractive. 
Using these values, the dimensionless interaction and hoping 
parameters are i = 37.8 and C = 0.23. 

The phase diagrams for repulsive (agp > 0) and attrac- 
tive (asp < 0) Bose-Fermi interactions are shown in Fig. l2] 
as a function of the boson and fermion chemical potentials. 
In both cases, a qualitatively similar structure emerges: The 
transition between the pure fermion and the phase where bo- 
son and fermions form a homogeneous mixture (mixed phase) 
is first order (thick solid red line). For the particular values of 
parameters chosen in Fig. plto describe the ^^^Yb-^Li mix- 
ture, we numerically find that the transition is weakly first 
order close to the origin (fip, iJi^s') = (0, 0), i.e., the chemi- 
cal potential values for which the slope of the free energy at 
p^ — Q changes sign are very close to the values of chemi- 
cal potentials at the transition. On the other hand, the transi- 
tion between the pure boson and the mixed phases (thin solid 
black line) is second order, i.e. continuous. This transition 
coincides with the locus of points where the system first de- 
velops a Fermi surface, i.e., fip — sign{gBp)p% — 0, and 
therefore fip > {jip < 0) for repulsive (attractive) interac- 
tions. In App.[B] we carry on an expansion for small fermion 
density which allows to establish the nature of the phase tran- 
sitions where the number of Fermi surfaces changes from zero 
to one. There, we argue that this transition is continuous be- 
cause of the scaling that the Fermi kinetic energy has with 
the fermion density in 3D, while it would be first order if 
the fermion would also move in strictly ID like the bosons. 



This result is also in agreement with the conclusion reached 
by Viverit et al. for Bose-Fermi mixtures in 3D |27|, where 
they show that phase separation between a mixed phase and a 
pure boson phase cannot be realized in 3D. 

Finally, the transitions between the vacuum, corresponding 
to zero density of both fermions and bosons, and either the 
pure boson or fermion phases (thin solid black lines in Fig. [2]) 
are continuous, as they correspond to the filling of a band 1 54 1 . 
Therefore, the first order line separating the pure fermion and 
mixed phases terminates at the origin (/i^, /i)/') = (0, 0) in a 
tricritical point (filled blue circle), where the first order transi- 
tion becomes second order. 

A first order transition in the phase diagram in chemical po- 
tential space implies that the system exhibits phase separation 
in density space, where, rather than fixing the chemical po- 
tentials jip and ^^^ , one fixes the boson ub and fermion np 
densities (see Fig.[3]l. We obtain therefore that, for finite inter- 
tube hopping t, phase separation is only possible between pure 
fermion and mixed phases (dot-dashed black lines in Fig. [3]l. 
In Sec. |IIIB] we will see that this fact is related to the dimen- 
sionality where fermions move and that the situation drasti- 
cally changes for strictly ID, e.g., when the hopping t for 
fermions is reduced to zero. Note also that, for attractive in- 
teractions, in contrast to a 3D Bose-Fermi mixtures in the ab- 
sence of the lattice l[T]|6], the system is found to exhibit phase 
separation rather than collapse. This result was also obtained 
in Ref. Ii25j . by treating the boson-boson interactions within 
the mean-field approximation. However, different from that 
work, the non-perturbative treatment of the boson interactions 
employed here, yields a first order transition between the pure 
fermion and mixed phases. 

The absence of a continuous phase transition between the 
pure fermion and any phase with a finite density of bosons 
can be qualitatively understood by making an analogy with the 
Landau theory of phase transitions and considering the series 
expansion of the free energy / in (|26| for small values of the 
boson density p^. In this limit, for fixed C,^^ ex gjj^, we 
have that 7 — {2Cp%)~^ ^ 1, and therefore the Bose gas 
is essentially fermionized and close to a Tonks gas. Thus, 
we can use the following asymptotic formula for the boson 
energy lEUEa: 



e(7) ~ CTG 

7— >-oo 



7 



+ 0(7-') 



where era — 7r^/3. Note that this expression implies that the 
boson contribution to the free energy grows as {p%Y. This 
yields the following series expansion for the free energy at 
small boson density: 



/ = /O + ~f2P% + UpI? + kiPlf + • • • , 

where the coefficients of the expansion are given by: 

/2 = sign{gBp)Ci ~ ii]P , 

h = -C_i/4 , 

/e = era - sign(5i3F)C_3/24 , 



(29) 
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FIG. 3 : (Color online) Phase diagrams in the density plane (ns , np ) ■ 
In order to obtain a more accurate estimate of the phase boundaries, 
we have used the parameters t — 0.023 and (" — 0.27. The corre- 
sponding phase diagrams in the chemical potential plane (not shown 
for brevity) display the same features and phase topology as the dia- 
grams shown Fig. [2] which are computed for the experimentally rel- 
evant a ^^^Yb-'^Li system, but for which the phase boundary in the 
density plane proved much harder to determine numerically. The top 
panel corresponds to a repulsive Bose-Fermi interactions, whereas in 
the bottom corresponds to attractive interactions. In both cases, the 
system can either be in a uniform mixed phase (lightly gray shaded 
region) or, by undergoing a first order transition, it can be in a phase 
separated state (dark gray shaded region), where a pure fermion and 
mixed phase coexist. The dot-dashed lines connect points on the first 
order boundary with the same values of the chemical potentials. 



and where C„ = 



Re 



y/jlp - e(k_L; 



> 



{\n\ is odd). In addition, note that C„ = if fip = 0. Thus, 
for jjp > 0, a continuous phase transition cannot take place 
because the coefficient /4 of the above expansion is always 
negative meaning that for /2 > the free energy must even- 
tually decrease away from the origin where p'g = before it 
can rise again (/g > is assumed, for stability). Thus, the 
free energy develops a local minimum for p'g ^ 0, which 
eventually can be tuned to be degenerate with the local min- 
imum at p^ = 0. It is worth comparing this situation with 
the result of the mean-field treatment of bosons interactions 
carried out by one of us in Ref. Il25l , where it was found that 
/4 = (2 — CC_i)/(4C), thus allowing both second and first 
order phase transitions to occur for /ij? > by tuning f2 — 
and /4 — 0. 

In addition, we can use the above expressions to understand 
the emergence of a tricritical point, which coiTesponds to the 
conditions f2 = = fi while /e > for stability. Since 
Cn = only for fip = 0, we can conclude that a stable tricrit- 
ical point can only exist at the origin of the chemical potential 
plane, i.e. for (/iFj/^s^) = (0,0). Close to the tricritical 
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FIG. 4: (Color online) Zero temperature phase diagrams for a repul- 
sive (aBF > 0, top panel) and attractive {ubf < 0, bottom panel) 
Bose-Fermi ^^*Yb-''Li mixture in ID, i.e., for zero hopping ampli- 
tude t = Q, and for zeta — 0.23. Differently from the finite t case 
(cf. Fig. [2]l, the transition between the pure boson and mixed phase 
becomes first order (thick red solid line), meeting with the other two 
first order transition lines (one between the pure fermion and mixed 
phase and the other between the pure fermion and the pure boson 
phases) at a triple point (filled violet square), where the three phases 
coexists. Bottom panel: For an attractive mixture there is neither a 
triple point nor a tricritical point, rather the first order line crosses 
the second order line at two critical end-points (filled [red] diamond) 
delimiting the region of phase separation between the vacuum and a 
mixed phase. 



point, the shape of the first order line can also obtained an- 
alytically using the conditions /2 > and /4 = ~2y f2fe- 
Note that for the choice of parameters done to describe the 
^^''Yb-^Li mixture, we find that, the true first order transi- 
tion obtained numerically stays very close to the one found 
analytically here. In addition, the transition is weakly first or- 
der because of the large values of the coefficient /g for which 

/2^0. 



B. Pure ID limit 

Next, we focus on the pure ID limit, i.e., the limit where 
the fermions, like the bosons, cannot hop between the tubes 
(i.e. t — 0). The phase diagram in chemical potential re- 
sulting from minimizing the free energy is shown in Fig. HI 
for both repulsive and attractive Bose-Fermi interactions. It 
can be seen that, in the pure ID limit, for both repulsive and 
attractive interspecies interactions, all transitions (except for 
the trivial ones from the vacuum phase) are discontinuous. In 
paiticular, the transition between the pure boson to the mixed 
phases, which was found to be continuous in the mixed dimen- 
sional system, becomes discontinuous as soon as the fermions 



are confined to ID. In App.JBJwe carry on an expansion for 
small fermion density that allows us to establish the nature 
of the phase transitions where the number of Fermi surfaces 
changes from zero to one. As shown there, the main difference 
between the mixed-dimensional case illustrated in the previ- 
ous section and the pure ID limit analyzed here can be traced 
down to the different scaling of the Fermi kinetic energy with 
the lineal fermion density p'^ in ID and 3D. In particular, 
whereas in ID the Fermi kinetic energy scales as {p'p}^, in 
3D it scales more slowly as (p^)^^"^- 

Furthermore, similarly to what was found in the previous 
section, the fermionization of bosons in ID also renders the 
transition between the pure fermion and mixed phases dis- 
continuous. These results are compatible with the previous 
results obtained by Das in Ref. IISTI . where was found that 
the transition between the pure boson and fermion phases is 
discontinuous thus leading to phase separation. 

The main difference between the repulsive and the attrac- 
tive cases is the way the transitions across which the density of 
fermions changes connect with the transition between the pure 
fermion and mixed phase. In particular, in the repulsive case 
(upper panel of Fig. |4]i, the three first order transition lines 
(thick red solid curves) between pure boson-pure fermion, 
pure fermion-mixed, and pure boson-mixed phases meet at a 
triple point (filled violet square symbol). At this triple point, 
the three phases coexist since the free energy exhibits three 
degenerate local minima (see Fig. E\. On the other hand, 
for the attractive case (see lower panel of Fig. HI, the triple 
point is absent. Instead, two critical end points (filled red dia- 
monds) appear. In the density phase diagram, the critical end 
points delimit a triangularly-shaped region (see bottom panel 
of Fig. l6]), where phase separation occurs between the vac- 
uum and mixed phases. Similarly to conclusion reached in 
Ref. ||25]| , this region can be regarded as a remnant of the col- 
lapse that occurs in the absence of a lattice in 3D Bose-Fermi 
mixtures with sufficiently large attractive interactions [1^6]. 

Let us finally remark that ID Bose-Fermi mixtures in the 
exactly solvable limit of equal masses (i.e. niB = rnp) and 
equal interactions strengths (g^ = gsp) have been ana- 
lyzed in Ref. |35|. By relying on a linear stability analysis, 
which requires that the compressibility matrix must be posi- 
tively defined for any strength of the interactions, the authors 
of Ref. 1*351 concluded that the system is always stable against 
demixing. However, in this work, by globally minimizing the 
system free energy, i.e., by looking for the global energy min- 
ima, we find that phase separation occurs in a rather broad re- 
gion of the density phase diagram and, in particular, it always 
occurs for small bosons and fermion densities in the pure ID 
limit. 



IV. CONCLUSIONS AND OUTLOOK 

To summarize, we have studied the equilibrium phase di- 
agram of a Bose-Fermi system in a mixed-dimensional ge- 
ometry of an array of coupled ID tubes, where bosons are 
strongly localized along the ID tubes, while fermions can 
hop between the tubes. Because we treat the boson interac- 
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FIG. 5: (Color online) Dimensionless free energy density / plotted 
versus the dimensionless bosonic density ub = p% for fixed di- 
mensionless hopping strength t = 0.023, fixed interaction parameter 
C, = 0.27, and a ID repulsive Bose-Fermi system, corresponding 
to the density phase diagram shown in the top panel of Fig. [6] The 
values of the chemical potentials are fixed at the three different first 
order transition lines near the triple point {p^ /Ers, pp/ Etb) — 
(35.3, 11.2), which describes the state where phase separation oc- 
curs between the three phases, pure boson, pure fermion and mixed 
phase: (1) ~ (36.3, 11.4) describe phase separation between pure 
fermion and mixed phase; (2) ~ (33.9, 10.8) phase separation be- 
tween pure fermion and pure boson; finally (3) ~ (38.1, 12.4) phase 
separation between mixed and pure boson. 
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FIG. 6: (Color online) Phase diagrams in the density plane (ns , np) 
for (^ — 0.27 and zero hopping strength t = (same remarks as 
for Fig. PI apply). Top panel: for repulsive interactions in the ID 
limit, phase separation occurs, not only between a pure fermion and a 
mixed phase, but also between a pure boson and pure fermion phase, 
as well as between a pure boson and a mixed phase. Bottom panel: In 
contrast, for the attractive case, there is no phase separation between 
pure phases, rather, it exists a region of phase separation between the 
vacuum and a mixed phase. 
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tion term non-perturbatively using the Bethe ansatz, we have 
found that the transition between the pure Fermi phase and 
the mixed phase is always first order This impHes that, for 
this system, phase separation can take place between a pure 
fermion and a mixed phase. However, the phase transition be- 
tween pure boson and mixed phases is found to be continuous 
for finite hoping amplitude of the fermions and phase separa- 
tion cannot take place between these two phases. This con- 
trasts with the results obtained by assuming that the bosons 
form a quasi-condensate and thus applying a standard mean- 
field treatment where the boson density is replaced with its 
expectation value ||25]| . In that case, the free energy the tran- 
sition between the pure fermion and a mixed phases can be 
continuous as well as discontinuous, depending on the system 
parameters. 

In the pure ID limit, we found that the all transitions are 
first order, except for the trivial ones between the vacuum state 
and the pure fermion or boson phases. The main differences 
introduced by the sign of the interaction are the existence of 
a triple point for repulsive interactions, and the appearance of 
two critical end points delimiting a first order line between the 
mixed and the vacuum phases. 

Thus, the main difference between the ID limit and the 
mixed dimensional system is the change in the character of 
the transition between the pure boson and mixed phases. We 
have argued (see App.JBJ) that this difference is a consequence 
of the different scaling of the Fermi gas kinetic energy contri- 
bution to the free energy function with the fermion density. 

Before commenting on the relevance of our findings for 
the experiments using ^^''Yb-^Li mixtures, it is interesting to 
note that, because we have applied the Bethe ansatz to the 
ID Bose Hamiltonian -ff]^ ( [T7] i, we can extend our analysis 
to to Fermi-Fermi system consisting of a light and a heavy 
atom. In fact, in the limit gBB -^ +00 where 7-^-00, the 
boson energy, 6(7) becomes identical to that of a free Fermi 
gas 1561 . By taking this limit of the Bethe-energy in Eq. (|26|, 
we found that the phase diagrams for both repulsive and at- 
tractive Fermi-Fermi interactions evaluated in the same way 
as before are qualitatively similar to those displayed in Fig. [2] 
for the Bose-Fermi mixture and that the differences are small 
and only quantitative |63|. 

The actual experimental systems are rendered inhomoge- 
neous by the existence of harmonic trapping. For mixed 
dimensionality systems, in general, and the ^^^Yb-^Li sys- 
tem ||441 -|46I loaded in an anisotropic optical lattice in par- 
ticular, we can rely on the local density approximation and 
deduce the main implications for experiments from the phase 
diagrams in the chemical potential space shown in Fig. l2] and 
in Fig. |4] For a given trap frequency and total atom numbers 
Np and Nb, it is possible to determine the range of values 
of the chemical potentials iip and fi^ of the phase diagrams 
in Figs. I2] and I4] being sampled by the trapped system. This 
range determines a region in the chemical phase diagram that 
contains the possible phases that will coexist in the trap. 

Lastly, we comment on the accuracy of the above mean- 
field approach that we have employed for the Bose-Fermi in- 
teraction term. Indeed, as any other mean-field theory, it ne- 
glects fluctuations. We expect fluctuations to be especially 



important close to a phase transition. Nevertheless, as shown 
later in Sec. [Ill] our calculations indicate that many of the 
phase transitions in the mixed dimensionality Bose-Fermi sys- 
tem (cf. Ill A I and in the pure ID limit are discontinuous. 



This means that, even very close to the transition point, fluc- 
tuations of the boson and fermion densities are typically sup- 
pressed and thus it can be expected that the mean-field the- 
ory to give a reliable picture. However, we also numerically 
observed that some of the transitions are only weakly first or- 
der In other cases, however, the transition was found to be 
continuous. Therefore a careful assessment of the effect of 
fluctuations will be required but will not be pursued here. 

Furthermore, within the above mean-field approach, we 
also have neglected the inter-tube couplings as well as the ef- 
fect of the bosons on the fermion properties, which are as a 
non-interacting Fermi gas. These are also concerns that de- 
serve to be investigated in the future work. Here, we have 
assumed that such effects are relatively weak and can only be- 
come important only for large values {ubfI and/or very low 
temperatures that may not be easily achievable under current 
experimental conditions. 

Beyond the assessment the accuracy of the present mean- 
field approach, another interesting direction is to apply the 
methods developed here to mixed dimensionality systems 
where the ID bosons interact via longer range interactions, 
such as dipolar gases or are tuned to the so-called the super- 
Tonks regime 120, 57-61 1 . 
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Appendix A: Mean-field approximation 

In this appendix, we provide the intermediate steps neces- 
sary to obtain the mean-field Hamiltonian, Eq. ( fTSj l. First of 
all, let us evaluate the expectation value of the fermion den- 
sity operator in a tube pptiix)- By substituting the expres- 
sion (|6]l into ( [TT| and using that, for a non-interacting Fermi 
gas, {.fl^f-^,) = (5k,k'?^k, with Tik being the Fermi-Dirac distri- 
bution function, we obtain: 



^i:- 



rik / dri|^R(ri)n0k,(ri)|2 . (Al) 



{PFIt{x)) 



Further, we can assume that the fermions occupy only the low- 
est Bloch band of the square lattice, thus the Bloch wave func- 
tion 0ki (r_L) can be developed in terms of Wannier functions 
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u'j^(rj^) localized at the R-th tube: 

0k, (r^) - ^ E e'^''-^^«R(r_L) , (A2) 

where N = NyN^. Therefore, Eq. ( |A1| | reads: 

(Pi.R(-)> = ^E"4 ^ j dv^^'^-^^'-^"-^ 



R'.R" 



X'u;^,(r_L)wR,„(r_L)|(^R(r_L) 



Nf 1 

IT TV 



j dv^\w^{v^)\^\^^{v^)\^^p°pA, (A3) 



where A'i? = ^j^ ny, p'p is the ID density of fermions ( fTS) , 
and 

^-A^/rfr^lV'R(r±)n0kJr^)P 

c.^y'dr^|«;R(r^)n^j,(r^)|2. (A4) 

In this derivation, we have used the fact that the boson Wan- 
nier orbital (^j^(rj^) is strongly localized around rj^ ~ R, 
thus we have neglected the contributions of the terms R' 
and R" different from R, because the corresponding wave- 
function overlaps are negligible. Note also that the constant 
A \AA) does not depend on the tube index R. Further, we 
can approximate ip-^{v±_) ~ e^'^'-^^^' /^^^/(a/tt^b) and 
wi^(rj_) ~ e-l""^-^l'/24/(y7^^^) with Ib < tp <i d, 
therefore the constant A is given by: 



A' 



n{£l+il) 



(A5) 



Taking into account that the mean-field averages of the bo- 
son and fermion densities are given by the expressions ( [T3] l, 
in order to obtain the final expression of the Hamiltonian, 
Eq. ( [T6| ), we need to deal with 

y^ / dx{pBRix))pFR{x) (A6) 

R -^ 

-PbY.J2 /rfr^l^R(r^)n'/'kJr^)P/^4 (A7) 

R k "^ 

^Ps^E/k/k- (A8) 



Appendix B: Series expansion for small Fermi density 

In this appendix we want to carry on an expansion for small 
fermion density so that to be able to establish the nature of the 
phase transitions where the number of Fermi surfaces changes 
from zero to one, such as the phase transition between a pure 
boson and a mixed phase. This is done in the same spirit to 
the expansion for small boson density conducted at the end of 
Sec, nil Al that has allowed us to establish that the transitions 



between pure fermion and any phase with a finite density of 
bosons can only be first order 

The starting point is the mean-field Hamiltonian derived in 
Eq. ([T6]l, with the difference that now, to simplify the anal- 
ysis, we will assume that the fermion dispersion appearing 
in Eq. ( fTSJ ) is quadratic, £(k) = k'^/2mp, where k = fc^; 
if fermions, like bosons, move strictly in ID, whereas k — 
{kx,ky,kz) if we instead consider the mixed-dimensional 
case where fermions are free to move in all three dimensions, 
while bosons still move strictly in ID. Note that assuming a 
quadratic isotropic dispersion for the fermions is expected to 
be a good approximation in the small fermion density regime 
that we are going to consider here, even if the true disper- 
sion in the mixed-dimensional case is anisotropic: In fact, 
when fermions start to mix with bosons, they must necessar- 
ily occupy the lowest band energy levels, for which the dis- 
persion can be approximated as quadratic. The anisotropy can 
be rescaled out when evaluating the fermion density p^ as a 
function of the Fermi wavevector kp, and yields to an overall 
prefactor relative to the isotropic result. 

Considering the mean-field Hamiltonian derived in 
Eq. ( [T6| ), the following step is to average over the fermion 
operator density in the fermion Hamiltonian ( fTS) , obtain- 
ing the Fermi-Dirac distribution function at zero temperature, 
(/k/k) = "-k = d{kF — k). The free energy density we ob- 
tain this way will have a different dependence on the Fermi 
density depending whether fermions move in ID or 3D. 

Let us start considering the case where fermions move 
in ID, so that the lineal Fermi density defined in Eq. ( [T5| ) 
is p^p = Np/NL — kp/n and the free energy potential 
g^ (i7™/)/LiV reads as: 



ip°pf - {^^F ~ gBFAp%) p%- 



6mp 



2m I 



-e(7) 



fi]fp%. (Bl) 



Note that now the free energy potential g — 
g{^^ , ^p, p'^,p^) depends on both the boson and fermion 
densities as well as on the chemical potentials. This means 
that minimizing g with respect to the fermion density, 

/(A^lf ,A*F,P^) = mmg{p]P,pp,p%,p''p) , (B2) 

Pf 

we retrieve the free energy potential /(/i)j^, pLp, p'g) consid- 
ered in Eq. ^T9\ , which global minimum with respect to the 
boson density gives the true thermodynamic grand-canonical 
free energy density ( |22] i. 

However, if we instead minimize g{p^ , fJ'F, p'b, p'p) with 
respect to the boson density. 



HPb ,Pf,Pp) ^myng{pg ,np,pg,pp), (B3) 

Pb 

SO that to eliminate p^ in favor of p^, p,p and p^, then 
the free energy potential h{p^ , pp, p^) thus obtained can be 
used to study the phase transitions where the number of Fermi 
surfaces changes from zero to one, by expanding for small val- 
ues of p'p. Because the dimensionless function e{j) coming 
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from the Bethe ansatz is determined numerically, this plan of 
deriving the free energy potential h{ii^ , fip, p'f) i^ not a triv- 
ial one. However, we can carry on this procedure in two op- 
posite limits, corresponding to high (7 = mBg]^/p^ ^ 1) 
and low (7 ^ 1) density, respectively: 




7< 1 
7> 1 



(B4) 



Note that the last limit results in a contribution to the free 
energy potential h of Tr^{p%)^ / 61113, which is identical to the 
kinetic energy of a free ID Fermi gas. This is the result of 
fermionization in the Tonks -Girardeau limit I20< . 

Evaluating ( |B3| l by solving dg/dp% = for p^, in both 
cases of high and low boson density, we obtain a series in 
terms of the Fermi density of cubic form: 



ho + h2p% + hi(p%f + he{p\ 



\3 

f) 



(B5) 



For high boson density the coefficients, for p^ > AgBpp'p, 
which requires p^ and gBP to have the same sign, are given 
by 



AgBF ID 

_^91f 
29Wb 



Pf 



Qmp 



(B6) 
(B7) 
(B8) 



In the low boson density limit, instead, we obtain (also for 
P^jP > AgBFPp) that the coefficients of the expansion ( |B5| l 
are now given by 



n 




Pf 



27r^2fi]f 

r^ ^^gBF^A^ 



6m I 



UV27Tp]P^^^ 



(B9) 
(BIO) 

(BID 



In this last case, note that, for m^ ^ niB, as we have as- 
sumed in the main text, the coefficient of the cubic term is ex- 
pected to be positive and large. In both cases, while /ig > 0, 
the coefficient of the square term is always negative, ^,4 < 0, 
implying that the transition cannot be continuous, rather it is 
first order. For small fermion density, the transition has place 
for /12 > and /i4 = —2\/h2hQ, implying that, for large val- 
ues of the coefficient hg, h2 —> and the transition is weakly 
first order, as we have indeed observed numerically in many 
cases. 

The p'^ expansion just carried on thus allows us to under- 
stand the nature of the transition between the pure Bose gas 



and mixed phases and deduce that it has to be first order in the 
strictly ID limit. We can merge this result with the one ob- 



tained in Sec. Ill A where we were expanding for small boson 



density (Tonks-Girardeau limit), and found that the transition 
between the pure Fermi gas and the mixed phases is also first 
order. We can thus conclude that in the ID limit the transitions 
between the pure boson/fermion phases and the mixed phase, 
must be all first order 

For comparison purposes, we can now consider the case of 
mixed dimensionality, i.e., where fermions move in 3D. Now, 
the 3D Fermi density is given by Np/^ = p'pd'^ = kp/6TT^, 
where the volume is i7 = NcPL, d is the spacing between the 
tubes (cf. FigfTb, and p^ the lineal fermion density. Thus the 
kinetic energy contribution to the free energy potential g = 
{H'^f)/LN now scales differently than in the ID limit: 



207T^m,F 



{pf - gBFAp%) p"p+ 
{P%f 



2mB 



e{l)~PBPB- (B12) 



As before, we eliminate the boson density p^ from the poten- 
tial g by making use of ( |B3[ ) and then we expand for small 
values of the Fermi density p^. In this case we obtain: 



(B13) 



h2p'F + hioMPpf^' 



where the coefficient ft,2 is given by the expressions ( |B6| l for 
the high boson density, while by ( |B9| l for the low boson den- 
sity limit. The coefficient /iio/3 is in both cases instead given 
by: 



^10/3 = 



^-4/3(g^2)5/3 

20TT^TnF 



(B14) 



We can thus see that, while the coefficient /i2 can change sign 
depending on the values of the system parameters, the coeffi- 
cient of the next order term, which scales Uke (p^)^^'^ rather 
than quadratically, is always positive /110/3 > 0. This means 
that the position of the closest minimum to p'^ = is entirely 
controlled by the sign of /12, that is, the transition between the 
pure boson and mixed phases must be continuous. By con- 



trast, as we have argued in Sec. Ill A the transition between 
the pure fermion and the mixed phase is independent on the 
dimensionality where fermions live and is always first order. 

The main conclusion of the simple exercise carried on in 
this appendix is that the different nature of the phase transi- 
tions where the number of Fermi surfaces changes from zero 
to one between the strictly ID case and the mixed-dimensional 
case is a consequence of the different scaling of the fermion 
kinetic energy in ID and in 3D. Note that this result is also ap- 
plicable to the transitions between the pure fermion and mixed 
phases, as bosons in ID at low density behave as a free Fermi 
gas by virtue of fermionization. 
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